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ABSTRACT 

A method of electric propulsion is explored where plasmoids such as spheromaks 
rnd field-reversed configurations (FRC) are formed and then allowed to expand down 
a diverging conducting shell. The plasmoids contain a toroidal electric current that 
provides both heating and a confining magnetic field. They are free to translate because 
there arc no externally-supplied magnetic fields that would restrict motion. Image 
currents in the diverging conducting shell keep the plasmoids from contacting the wall. 
Because these currents translate relative to the wall, losses due to magnetic flux diffusion 
into the v/all arc minimized. During the expansion of the plasma in the diverging cone, 
both the inductive and thermal plasma energy are converted to directed kinetic energy 
producing thrust. Specific impulses can be ir. the 4000 to 20000 sec range with thrusts 
from 0.1 to 1000 Newtons, depending on available power. 


BACKGROUND 


While chemical rockets are capable of veiy high thrust, consumption of propellant 
is excessive for many space missions such as long-term satellite station keeping and 
planetary missions. Their nozzle exhaust velocities are limited by heats of reaction t,o 
~ 5000 m/sec, limiting specific impulse to I, ~ 500 sec. Nuclear thermal rockets 
using hydrogen propellant can approach 950 sec, still too low for many missions. 

For many space missions, specific impulses in the 2000 to 10000 sec range are 
desirable in order to limit propellant consumption. However, me power in the exhaust 
stream 


P e = ~ rhul , 

must be supplied by means other than chemical. Here rh is the propellant mass flow 
in kg, sec. The only practical form is electrical power and an electric generator, either 
solar * r nuclear (fission or fusion), must be carried on board. (Batteries, fuel cells, 
and other chemical energy systems are not. suitable because their Am! supply or energy 
storage mass is comparable to, and subject to the same limitations a.s, chemical rockets.) 
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When one trades off power plant mass against propellant mass, it becomes clear that 
every mission has an optimum specific impulse above which power plant mass becomes 
excessive and below which propellant mass dominates. 

In terms of piopulsion, space missions can be characteris'd by a certain payload 
mass m p carried through a velocity increment AV. Typical AVs range from S km,/ sec 
to go from low earth to geosynchronous earth orbits (LEO - GEO), to 14 kn /sec *■:> go 
from LEO to Mars orbit and return, to 110 km/ sec to go from LEO to Saluki orbit 
and return (ref 1). Pay! ad masses can range from 10s of kg lor a small experimental 
satellite to hundreds of tonnes for a. manned Mars mission. 

Existing electrically- powered thrusters can supply high specific impulse at rather 
low thrust. There are three kinds: electrothermal (resistojets and arejets), electrostatic 
(ion thrusters), and electromagnetic imagnetoplasmadynair’ic (MPD) thrusters]. 

Resistojets use resistance wire to heat a flowing gas. Arejets do the same with 
arc. discharges across *he grs Ion thrusters use perforated accelerator plates at high 
potential difference to accelerate an ionized gas. MPD thrusters are similar to arejets 
except they also exploit the J X B force in the arc to blow the plasma through a nozzle. 

Only the MPD thrus*er is capable of both high thrust and high specific impulse. 
However, these thrusters have rather short lifetimes due to electrode erosion (ref. 2). 
The plasmoid propulsion thruster discussed below may fill the need for high-power levels 
with long life and high specific impulse. 


DESCRIPTION OF THE CONCEPT 

The basic idea is shown in Fig. 1. Compact tori (CT) axe formed inside an 
electrically conducting chamber. A compact torus is a toroidal-shaped ionized gas con- 
taining an internal eiectric current traveling mainly along the minor axis of the torus. 
Formation can be either by a coaxial Marshall gun or by a rapidly pulsed induction 
coil. Because the internal current can twist helically at locations off the minor axis, the 
configuration could take two forms: (1) a spheromak. which has near-equal magnetic 
field strength at the plasma edge in directions parallel and perpendicular to minor axis, 
or (2) a field- reversed configuration (FftC) where the current has no helical twist and 
all of the magnetic field from that current is perpendicular to the minor axis. Both con- 
figurations are being explored for nuclear fusion applications. The difference between 
the two is that the FRO may be capable of higher plasma pressure but appears to have 
poorer confinement of the plasma thermal energy' than the spheromak. 

Once formed, the compact torus or plasmoid is pushed into a diverging conducting 
shell where it spontaneously accelerates down the divergence so as to reduce its internal 
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Fig. 1. Plasmoid propulsion concent. 

pressure. Image currents in the shell keep the plasma off the wall anti provide the 
physical push that gives the thrust. 

This is a pulsed system and thrust variation is accomplished by simply varying 
the rate at which these piasmoids axe formed. We believe that realistic pulse rates can 
vary from 100 Hz on down. 


ANALYSIS 

In this section, we examined the dynarmes and energetics of plasmoic. ..Apansion 
down the conducting shell in order to arrive at preliminary estimates of specific impulse, 
thrust, and losses. 

The magnetic energy in the plasma is 


Em — - LI’ 




where I is the internal current in the plasma. L is the plasma inductance and is 
very nearly /x 0 R lor aspect ratios R/a 1.5, typical of spheromaks. Here is the 


L r 
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magnetic permeability of free space; and a and R are the plasma minor and major radii, 
respectively. The thermal energy, assuming electron and ion temperatures are equal, is 


E t — 3n T V p , 


where n and T are the the ion density and temperature, respectively, and V p the plasma 
volume (— 2n 2 Ra 2 ). 

During an adiabatic expansion, total rmgnetic flux is conserved. That is, 


B p aR = constant , 


where B p is the magnetic field due to the plasma current I: 


B p = 


2 rc a 


Assuming geometric similarity during expansion (R ~ a), it can be shown that 


la -- constant , 


Therefore, as the plasma expands, the hitern<d total current drops linearly with size. 

Over an expansion from ai to a 2 , it can be shown fr m the above that the magnetic 
energy scales like 


Em2 _ «1 

Enal 


A five-fold increase in size will therefore com -t 8U% of the magnetic energy to kinetic 
energy of the partichs by interaction of the plasma current ana the image current in 
the inclined shell. 


The thermal energy in the plasma is also converted to kinetic 
expansion. With an adiabatic expansion 


energy during 
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Since, for hydrogen, 7 = 5/3, the thermal energy scales like 


E(2 _ ai 
Etl Q>2 


A five-fold linear expansion will 
plasma and, like *ne expansion in 


therefore recover 96% of the thermal energy in the 
a regular nozzle, convert it to kinetic energy. 


so of interest is the behavior of (3 during expansion. 3 is the ratio of plasma 
PreSSUre "" " 8 meaSUr ' efficiency of the 



It can be shown from the above equations that 
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where M is the total plasma mass: 

M = (1.67 • IQ -27 kg) n V p W 

Kere „ is the average ion density and W is the ion atomic weight. 

The primary issues which drive the ci 7 : n „ „r *u 
loss from the hot plasma (typically over 10000K) 'and ' 7” ‘j ‘' m ' f ° r rad, “ t!on 
relative to the acceleration time of the plasmofrh ’ ' time ” b ^ h 

Radiation calculations are verv comnW i , 

here. Hopkins and Rawls (r**f 3) have nr A a T* ^ Cruc ^ e est iniates can be made 
loss mechanisms. Figure 2 show ^ t0 the Vario - radiation 

total radiation loss, in IV/m^Ts b" 1 °" * fcr eleme "‘ s - 

Pr = Ip n e m , 

where n e and n* are the electron and ion density of the l 

characteristic radiation time is ty f * h pIasma ’ respectively. The 


the problem right^quires^drtSd't- 01 ° f ioss lhrou Sh radiation. To do 

assessing fcjfc* “"T- h “ * 

from the other elements. While it would be nir t ro ™ h >' dro S en is well under that 
raise thrust and lower snedfic Iml.T t * USC thc ’ ore elements to 

option and hydrogen is the element of droice'!"'’ °“ l lha ‘ '‘'“ ,li ° n loss P reclud ' s ‘his 

f or a gi™ pwa^r^r ** density 

.sat, on occurs. Figure 3 shows ionization of hydrogen 1 ThUSt ““l? part,al Ion - 
and density taken from ref. 4). At 1.5 eV and 10'®™“ f foncfron of temperature 

hydrogen is about 50% ionized Major concerns 1 tZ ^ *" ' “”P>« • 

ionization to occur Mid the interaction of the ions with thrre'''' 0 '"^ ** “f leTd ° f 
issues will be addressed in future work. “ remaining neutrals. These 

time U h U k tlrl“m ^-istive decay time, which is the 

resistance. If this is shor7como^d d,SS ' P, “ e dl " ‘° lhe “‘“"a 1 pWa 

raise the plasma temperature and dissipate' t le't'me'nT The’ “ ““ 

to the point of disruption. It is not clear if tins is a 
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Fig. 2. Total radiation loss parameter as a function of plasma electron temperature. Corona! 
equilibrium is assumed (from ref. 3). 

still expand down the nozzle converting its thermal energy to kinetic. And since, for a 
given expansion, a greater fraction of thermal energy is converted than magnetic, thrust 
efficiency could actually be higher. We will also explore this issue in follow-on work. 

The full expression for plasma inductance is (ref. 5) 


= A 1 qR 


2+ l i 




Fig. 3. Production of electrons in a hydrogen plasma as a function of particle density and 
electron temperature (from ref. 4). 


where U is the plasma internal inductance. For uniform current profiles, which is likely 
here (due to turbulence), l{ = 1/2. With this and an aspect ratio R/a = 1.5, the bracket 
term is very nearly l.C and so L ~ /x 0 R . 

The plasma resistivity is 

■q = 5.38 • 10~ 5 InAZ T~ ir> , 

f 

l 

: where Z is the atomic number of the gas and In A is a constant that is typically about 

16. The plasma total resistance, assuming flat temperature profiles and no resistivity 
enhancement due to turbulence or helical pitch of the current, is given by 


n = *7 


2tt R 
i r a 2 


| 

! The time for the current to decay to 1/e of its initial value is L/fl. 
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The characteristic times above must be compared to the time to accelerate the 

nl AJimnirl Hnnm of 'TV. K»r.U J.'.l _ d i i . « 

r ***** xv xiaxLiii xcMxiiu iuii6bit energy tue nozzie lengtn 

L n should be at least five times the outer radius of the plasmoid. Ass umin g constant 
acceleration, the time for acceleration is 


2Ln 

«e 


The above expressions are sufficient to perform an initial assessment of pulsed plasmoid 
propulsion. This is done L°low. 


SCOPING CALCULATIONS 


Using the above plasma equations along with some expressions for missions, a 
short code was written and exercised to determine operating parameter spaces. Two 
example cases are reported here. If time permits in the future, sensitivity studies could 
also be performed with this code. 

Code output for a small thruster is shown in Table 1. Inputs 1 through 7, 10, 
and 13 define the plasma while Inputs 8, 9, and 12 refer to the mission. The power 
plant specific weight of 2 kg/kW is not unreasonable for an advanced nuclear electric 
plant or an advanced photovoltaic solar array (ref. 0). The radiation parameter of 
7-10 38 W-m 3 is taken from Fig. 2 with hydrogen at 1.5 eV. Note that this temperature 
is used to establish the density at a given 0 . In actual fact, both temperature and 0 will 
float depending on the dynamics of the plasma power balance. That analysis is very 
complex and well beyond the scope of this effort. 

With plasma current, 0 , temperature, field, and aspect ratio dictated, plasma 
dimensions, density, total mass and energy are determined from the equations discussed 
above. Also determined are the radiation loss and resistive decay times, the acceleration 
times, specific impulse and thrust. Using the mission input parameters, one can readily 
determine power plant and propellant mass, and a lower bound on acceleration time 
to the specified AV (it is an underestimate for guidance only because just the payload 
mass is considered). 

The example in Table 1 is for a thruster that could be used to transfer the 1000 kg 
satellite from LEO to GEO in about 200 days. This is typical of electric propulsion 
missions and care is taken to protect delicate parts during the long spiral transit through 
the Van Allen radiation belts. With a specific impulse of 8000 sec and a thrust of 0.5 N, 
the total mass of propellant used is only 11.4% of the payload mass. Similarly, with an 

exhaust power of 20 kW, the power plant mass, here likely to be solar cells, is only 5% 
of the payload. 
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TABLE 1 

SMALL HYDROGEN THRUSTER FOR ORBIT RAISING 


Input List 


1 Plasma current, amps 

0.10E+06 

2 Poloidal beta 

0.20 

3 Plasma temperature, eV 

1.50 

4 Poloidal field, tesla 

0.10 

5 Aspect ratio, R/a 

1.50 

6 Ion atomic number 

1.0 

7 Ion atomic weight 

1.00 

8 Mission delta- V, m/sec 

9000 

9 Mission payload, kg 

1000 

10 Radiation parameter, W-m 3 

0.7CE-35 

11 Electron density fraction 

0.30 

12 Power pKnt specific weight, kg/kW 

2.0 

13 Pulse frequency, Hz 

10.0 

Output 


Specific impulse, sec 

8027 

Thrust, N 

0.515 

Ratio of radiation to accel time 

23.3 

Plasma acceleration time, sec 

0.636E-04 

Radiation time, sec 

0.149E-02 

Propel! ant -to-pay load mass ratio 

0.114 

Exhaust power, W 

0.203E+05 

Power plant mass, kg 

0.507E+02 

Ratio of power plant mass to payload 

0.051 

Total plasma energy, J 

0.203E+04 

Plasma thermal energy, J 

0.141E+03 

Plasma magnetic energy, J 

0.188E+04 

Plasma mass, kg 

0.655E-06 

Exhaust mass flow, kg/sec 

0.655E-05 

Ion density, m~ 3 

0.166E+22 

Electron density, m~ 3 

0.497E+21 

Plasma volume, m 3 

0.236 

Major radius, m 

0.30 

Plasma radius, m 

0.20 

Plasma o.d., m 

1.0 

Accelerator cone length, m 

2.50 

Plasma current resistive decay time, sec 

0.536E-04 

Spacecraft acceleration time, days 

202 
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Considerable iterations were required to find a configuration that did not radiate 
excessively. It turned out that magnetic fields must be low, 0.1 T in this case, resulting 
in fairly large plasmas with low density. This is no surprise since radiation power density 
goes like n 2 . In this case, the radiation time is 23 times longer than the acceleration 
time, providing a comfortable margin to compensate for the crudeness of the calculation. 
The plasma L/ft time of 53.6 /xsec, however, is comparable to the acceleration time of 
63.8 ^isec. This suggests that some magnetic energy will convert to thermal, and then 
to kinetic, during the expansio . Clearly, detailed analysis of all this is needed, but it 
is very complex. 

The overall plasma diameter of 1.0 m is quite large. However, although formation 
coils would have similar dimensions, their fields should be on the order of the magnetic 
field (0.1 T) and therefore would have fairly thin windings. The accelerator cone is 
2.5 m long with an maximum diameter of perhaps 5 m. The forces on this cone are 
very small and so it could perhaps be made of very thin aluminum. Such a cone made 
of 10-mil thick aluminum would weigh about 60 kg (6% of payload) and could double 
as a heat radiator (more on this later). 

A larger thruster, suitable for a manned Mars mission, is shown in Table 2. To 
increase thrust, plasma current is raised from 100 kA to 1.0 MA. Since thrust goes like 
I 3 , it increases from 0.5 to 515 N. Using 20 km/sec as a conservative AV for a Mars 
mission, a 200 tonne payload requires roughly 90 days worth of acceleration. Actually, 
it is longer because the propellant and power plant, each 25% as massive as the payload, 
must be added in as must other components discussed below. 

A somewhat fanciful outline of such a spacecraft is shown in Fig. 4. As seen in 
Table 2, the plasma has a 10 m overall diameter and an accelerator length of 25 m. 
With a 5:1 expansion, the cone outer diameter is 50 m. This is acceptable only if it can 
serve multiple purposes. It turns out that it is just about the right size for the reject 
heat radiator for the nuclear reactor. The reactor power conversion is assumed to be a 
helium closed cycle gas turbine with a mean heat rejection temperature of 650K (1000K 
cooling to 300K). Assuming an emissivity of 0.8 and a 30% power conversion efficiency, 
then the 58 MW(th) reject heat cam be radiated at roughly 8 kW/m 2 . The total radiator 
area must then be around 7000 m 2 . The above cone area is about 9000 m 2 , just a little 
more than needed, which is fine because it allows for added heat input from plasma 
radiation (which is then radiated directly to space). If the specific weight of the nozzle 
can be held to 5 kg/m 2 , its mass would be about 45,000 kg, which is an acceptable 22% 
of the payload mass. 

There is no reason why the accelerator nozzle must be conical. It is only necessary 
that the plasma expansion be smooth. If instead it is made parabolic, then it may also 
be able to serve as a high-gain antenna. A swing-boom signal collector can be placed 
at the focus when the thruster is off. During the acceleration period from earth, the 
antenna ; s pointing properly to receive terrestrial signals. 
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TABLE 2 

LARGE HYDROGEN THRUSTER FOR MARS MISSION 


Input List 


1 Plasma current, amps 

0.10E+97 

2 Poloidal beta 

0.20 

3 Plasma temperature, eV 

i.50 

4 Poloidal field, tesla 

0.10 

5 Aspect ratio, R/a 

1.50 

6 Ion atomic number 

1.0 

7 Ion atomic weight 

1.0 

8 Mission delta-V, m/sec 

20,000 

9 Mission payload, kg 

200,000 

10 Radiation parameter, W-m 3 

0.70E-35 

11 Electron density fraction 

0.30 

12 Power plant specific weight, kg/kW 

2.0 

13 Pulse frequency, Hz 

10.0 

Output 


Specific impulse, sec 

8027 

Thrust, N 

515 

Ratio of radiation to accel time 

2.34 

Plasma acceleration time, sec 

0.636E-03 

Radiation time, sec 

0.149E-02 

Propellant- to-pay load mass ratio 

0.254 

Exhaust power, W 

0.203E+08 

Power plant mass, kg 

0.507E+05 

Ratio of power plant mass to payload 

0.253 

Total plasma energy, J 

0.203E+07 

Plasma thermal energy, J 

0.141E+06 

Plasma magnetic energy, J 

0.188E+07 

Plasma mass, kg 

0.655E-03 

Exhaust mass flow, kg/sec 

0.655E-02 

Ion density, m“ 3 

0.166E+22 

Electron density, m -3 

0.497E+21 

Plasma volume, m 3 

236 

Major radius, m 

3.0 

Plasma radius, m 

2.0 

Plasma o.d., m 

10.0 

Accelerator cone length, m 

25.00 

Plasma current resistive decay time, sec 

0.536E-02 

Spacecraft acceleration time, days 

89.9 
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Fjg. 4. Schematic of space station for manned Mars mission. Note that the large rocket nozzle 
serves multiple purposes. 


DISCUSSION 


The major attraction to this thruster is the potential for high thrust and high 
specific impulse and, if electrodes are not needed for formation, long lifetime. One 
would also expect the energy efficiency, which is the ratio of exhaust power to power 
supplied, to be quite high. In these respects, plasmoid propulsion appears to fill a gap 
in the arsenal of thrusters now available. 

The main problems stem from radiation loss and the resulting large plasma 
dimensions needed to keep it in check and the partial ionization of the gas. The large 
size is compensated for by the low magnetic fields and low forces on the accelerator 
nozzle, permitting it to be made of thin- walled material. If the nozzle can in fact double 
as a heat radiator, then the large size is not a handicap. 

There has been no discussion thus far of the power conditioning required. This is 
a pulsed system and high voltages with fast risetimes are required to form the plasma. 
In the case of the Mars thruster, total energies around 2 MJ must be supplied. If 
capacitors were used in this case, their mass would be roughly 40,000 kg, about 20% of 
the payload mass. This may actually be tolerable, although other pulsed energy sources 
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such as inductors should be explored as well. If the pulse frequency can be increased, 
the energy per pulse will drop, lowering not only the capacitor mass but the mass of the 
large accelerator cone. The 10 Hz chosen for the examples seemed reasonable from the 
standpoint of circuit recharging and chamber clearing. It can probably be much higher. 

If one adds up all the major hardware for the Mars vehicle above, including the 
capacitors, one gets a mass roughly equal to the payload mass; i,e., the payload fraction 
is 50% of the total initial mass, which would be very appealing for a Mars mission. 


FUTURE WORK 


The next step is to perform time-dependent plasma analyses and to scope out 
the design of the power conditioning system and plasma startup system. All of this is 
needed before an accurate estimate cam be made of thrust efficiency, which is a very 
important number because it determines the size of the power plant needed. 
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